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ABSTRACT 



We investigate the variation of galaxy clustering with luminosity using 
the recently completed SSRS2 sample. Clustering measurements based on 
the two-point correlation function and the variance of counts in cells reveal 
a strong dependence of clustering on luminosity for galaxies brighter than 
L*, while no significant variation is detected for fainter galaxies. We derive a 
relative bias versus magnitude relation which can be compared with theoretical 
predictions. Existing models of galaxy formation cannot adequately reproduce 
the simultaneous steep rise of biasing at high luminosities and the plateau at 
the low-luminosity end. Improved modeling of halo-galaxy relation, and larger 
samples including low luminosities galaxies are required to draw more definitive 
conclusions. 



Subject headings: cosmology, galaxies: clustering, large-scale structure of the 
universe 
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1. Introduction 

Redshift surveys of galaxies have provided over the years a wealth of information that 
has greatly contributed to our understanding of the way galaxies are distributed in space. 
However, these observational results on galaxy clustering are not easily used to constrain 
theoretical models because of our poor understanding of the relation between galaxy and 
mass distributions. 

Although until recently the main emphasis of redshift surveys has been on the study 
of large-scale structure, several attempts have been made to understand the dependence of 
galaxy clustering on their internal properties such as luminosity, morphology, and surface 
brightness. An understanding of the relative bias between different galaxy populations is 
essential for constraining models of galaxy formation. Most of these studies were based on 
the first generation of redshift surveys hke CfAl (Huchra et al. 1983) and SSRS (Southern 
Sky Redshift Survey) (da Costa et al. 1988) which were not suitable for this kind of analysis 
because of the small volume surveyed and of the relative small number of galaxies in the 
samples. 

Not surprisingly the results have so far been rather controversial, with some authors 
claiming evidence of luminosity segregation in the data (see for instance Hamilton 1988, 
Davis et al. 1988, Borner, Mo and Zhou 1989, Domingucz-Tcnreiro & Martinez 1989, 
Borner, Deng and Xia 1989, Maurogordato & Lachieze-Rey 1991, Maurogordato, Schaeffer 
and da Costa 1992) while others arguing against it (as PhiUipps & Shanks 1987, Thuan 
et al. 1987, Alimi et al. 1988, Eder et al. 1989, Haynes & GiovaneUi 1989, BingeUi et al. 
1990, Thuan et al. 1991, Hasegawa & Umemura 1993). While the dependence of clustering 
on morphology is well estabhshed (Davis & Geller 1976, Santiago & da Costa 1990, Borner 
& Mo 1990, Mo et al. 1992, lovino et al. 1993, Dominguez-Tenreiro et al. 1994), the 
dependence on luminosity is still a matter of debate. 
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Recently, a new generation of redshift surveys has been completed calling for a new 
examination of the problem. Recent estimates of the galaxy power spectrum from the CfA2 
(Park et al. 1994), and of the correlation function for Stromlo-APM (Loveday et al. 1995) 
redshift surveys seem to show that some luminosity segregation does exist, but a consensus 
has still not been reached over the range of magnitudes at which the effect occurs. In this 
paper we will analyze in detail the dependence on luminosity of the correlation function 
and of the variance of counts in cells using the recently completed SSRS2. The combination 
of dense sampling, large sky coverage and a large number of galaxies allows us to define 
various subsamples to investigate separately different effects that may have a bearing on 
our conclusions. In section 2 we describe the sample and the corrections applied to the 
data. In section 3 we discuss the different statistics applied in our analysis, investigate the 
dependence of clustering on luminosity and compare our results to those of previous work. 
In section 4 we discuss the implications of our results on models of galaxy formation. A 
brief summary of our main conclusions is presented in section 5. 

2. The data 

In the analysis presented below we use the SSRS2 (da Costa et al. 1994), a 
magnitude-hmited redshift survey centered on the south galactic pole. It includes ~ 3600 
galaxies and covers 1.13 sr, delimitated by —40° < 6 < —2.5" and b < —40°, complete up to 
limiting magnitude 1713(0) — 15.5. Radial velocities have been corrected to the Local Group 
rest frame (Yahil et al. 1977). 

Since all galaxies in the SSRS2 have been assigned morphological types following 
the ESO-Uppsala classification system (T), we have adopted a morphological-dependent 
K-correction, K{z,T), to compute the absolute magnitude of a galaxy according to the 
expression 
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M ^m-25-5logDL- K{z,T) (1) 
where Dl is the luminosity distance for standard Friedman models: 

Dl = ^^[1 - ?o + qoz + (go - i)\/'2qQZ + 1] (2) 

-no Qo 

In the following analysis we have adopted Hq = 100, qo = 0.5. We have divided the 
sample in four morphological classes and applied for each class the K-correction following 
Efstathiou, EUis & Peterson (1988): 



E/SO -5 < T < 4.Uz 

Sa, Sb, Sbc 1 < T < 4 2.90z 

Sc, Scd 5 < T < 6 K = 2.25z 

Sdm 7 < T K= 1.59z 

Volume-limited samples were defined adopting a "minimum distance" criterium by 
taking galaxies within a depth corresponding to the smallest among those calculated 
for the different morphological classes, which corresponds in our case to the depth of the 
early-type galaxies. Luminosity distances were then converted into comoving distances. 

3. Analysis 
3.1. Statistics 



We examine the clustering properties of galaxies using the two-point correlation 
function ^(s) and the J3 integral, computed by counting the number of galaxy pairs 
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in shells, and the variance of galaxy counts in excess of Poisson derived from count 
probabilities in randomly placed shperes. 

The two-point correlation function in redshift-space was computed using the following 
estimator (Hamilton, 1993): 



where DD is the number of the galaxy-galaxy pairs, DR is the number of galaxy-random 
point pairs and RR is the number of pairs in the random catalog with a separation in the 
interval s and s ^ ds. This estimator is less sensitive to uncertainties of the mean density, 
which is a second order effect in contrast to the previously defined estimator (Davis & 
Peebles 1983). Since we only consider samples limited both in distance and in absolute 
magnitude we adopt uniform weighting in the pair-counts defined above. From these pair 
counts in shells, we have also computed the so-called J3 integral (Peebles 1980) 



which represents (up to a factor of 47r) the mean number of excess galaxies around each 
galaxy within a distance s. We use this statistic in section 4 to derive the relative bias 
between different luminosity classes. 

Another way to characterize the clustering properties of galaxies consists in calculating 
the probability P(A^, /) of finding N galaxies in spherical cells of radii / placed at random. 
In this procedure, the spheres crossing the edges of the samples are rejected. The number 
of trials is chosen in order to have small statistical errors at scales larger than IhT^ Mpc 
(see the discussion in Appendix B of Maurogordato, Schaeffer, da Costa, 1992). This leads 
in our case to 10*^ to 10^ trials, from the shallower to the deeper sample. We have derived 
the variance of these counts in cells which can be expressed as 




(3) 




(4) 



< (AT- Ar)2 >^ iV(l + Ar,), 



(5) 
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where n is the mean density of the sample, N = nV is the mean number of objects, and 
Nc = n'^V'^a'^ is the mean number of objects in excess of random inside a sphere of volume 
V. The variance can be related to the volume-averaged correlation function 

a\s)^as) = ^^j^as)d's,dh2 (6) 
In the power law approximation, ^(s) = (so/s)"^, this can be expressed as 

^'<-'> = (l-oV3)(l!oV4)(l-oV6)«-'> 
In this case ^3(5)/^ and cr^(s) are related by a multiplicative factor depending weakly on 7 
in the range of interest (0.8 for 7 = 1, 0.75 for 7 = 2). Below we will compare the results 
obtained from both methods to test our conclusions. 



3.2. Volume and Luminosity limited Samples 

From the SSRS2 we have extracted nine volume and luminosity limited subsamples 
with their absolute magnitude limit separated by 0.5™, ranging in depth from 30 to 168 /i^ 
Mpc . The characteristics of these samples are given in Table 1 where we list: in column 
(1) sample identification; in column (2) the depth as defined in the previous section; in 
column (3) the corresponding faint absolute magnitude limit of the sample; in column (4) 
the number of galaxies in the sample; in column (5) the mean density of galaxies. 

EDITOR: PLACE TABLE |1] HERE. 

For each of these samples we have computed ^(s), ^3(5) and cr^(s), and derived 
power-law fits to S,{s) in the range of separations where it is a good approximation (this 
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range is not the same for all the samples). The correlation parameters derived from these 
fits are listed in columns (6) and (7) of Table 1 where we give the correlation length sq 
and slope 7 for a two-parameter fit in the separation range range given in column (9). 
For comparison we list in column (8) the correlation length Sq obtained assuming a fixed 
slope 7 = 1.7. In column (10) we give the value of the variance in 8h^^ Mpc spheres, o"| as 
derived from the analysis of counts in cells. Finally, in column (11), we list the ratio L/E 
of late type over early type galaxies. 

EDITOR: PLACE FIGURE |l| HERE. 

Figure 1 shows Ci^)} Jsi^) ^iid cr^(s) for the samples given in Table 1. The error bars 
correspond to the rms value as measured for 15 bootstrap re-samplings of each sample. In 
the case of the less dense samples, errors on the various statistics on small scales can be 
very large. Therefore, here and in the following graphs, error bars for these samples are 
shown as dotted and dashed lines in order to avoid over-crowding of the plots. While we do 
not find any systematic variation of the slope with depth (table 1), the correlation length 
varies significantly over the range of depths considered especially for the deepest samples 
containing the brightest galaxies. The observed variation is consistent to that shown in 
Davis et al. (1988) but now probing much larger depths for a given absolute magnitude 
limit. The same trend for an apparent increase in the clustering strength with depth is seen 
in the number of neighbors Jsis) and in the variance computed from counts-in-cells. This 
variation is quantified by (t| listed in table 1. 

A key question is whether the observed variation reflects the intrinsic nature of the 
galaxy distribution (e.g. Baryshev et al. 1995 and references therein), simply sampling 
variations or genuine luminosity segregation, since in deeper samples the fraction of 
intrinsically luminous galaxies increases. The discussion is particularly problematic since 
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comparison of the amplitude of the correlation function as determined from different 
volumes may be affected by sampling variations which can arise because different structures 
are probed for different volume-limited samples. This could induce fluctuations in the mean 
density of the sample or variations in redshift distortions as the sample probes different 
virialized systems and samples different regions of large-scale peculiar velocity field. Still 
another effect that could be introduced by redshift distortions would be if more luminous 
galaxies were found preferentially within virialized systems, as a product of mergers for 
instance. In this case, redshift distortions could also affect the comparison of samples even 
within the same volume. 

Ideally, to prove the existence of luminosity segregation one should compare the 
clustering properties of different luminosity classes within the same volume and in real 
space. Recently, this has become partially possible thanks to the large number of galaxies 
in complete dense surveys like the SSRS2. The large number of galaxies allows the sample 
to be divided in different ways to discriminate between the various effects discussed 
above. This can be achieved by comparing the clustering properties of galaxies of different 
luminosity classes within the same volume and the clustering of galaxies in the same 
luminosity class but in different volumes. Although dividing the sample in this fashion 
reduces the number of galaxies in each subsample, thus increasing the uncertainties of the 
various statistics, it makes it possible to distinguish the contribution that different effects 
may have on the ^(s) derived from volume-limited samples like those used in figure 1. 

EDITOR: PLACE FIGURE | HERE. 

In figure 2 we show $,{s) for three different volumes {D = 38h~^ Mpc, D = 60h^^ 
Mpc, D = 91h^^Mpc). For each volume we have compared the clustering properties of 
galaxies in contiguous magnitude bins 1"* wide. Unfortunately, the sample is not large 
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enough to allow this comparison to be made over a wider range of luminosity classes since 
we are constrained by the apparent magnitude limit of the sample and by the number of 
galaxies in each class. In our case we are limited to only two adjacent magnitude bins given 
in Table 2. The depth was chosen to correspond to the maximum distance at which the 
faintest galaxies considered are visible. Dividing the sample in this way allows us to test for 
luminosity segregation in a way that is insensitive to sampling variations. 

EDITOR: PLACE TABLE | HERE. 

We note from figure 2 that faint galaxies do not exhibit any significant relative 
bias. While there is some slight indication of segregation for galaxies in the range 
—20.5 < M < —19.5 relative to galaxies in the magnitude interval —19.5 < M < —18.5 
(panel b), it is only in the largest volume considered that we find a significant effect. This 
result implies that galaxies brighter than —20.5 are significantly more clustered than those 
in the range —20.5 < M < —19.5. Based on these results we infer that the onset of this 
segregation occurs roughly at M = -19.5, close to L*. 

EDITOR: PLACE FIGURE | HERE. 

In order to examine the importance of sampling variations we show in figure 3 
measurements of S,{s) for galaxies in the same luminosity class in three different volumes 
of increasing depth. The characteristics of these subsamples are also summarized in Table 
2. From the plot we find that the clustering of faint galaxies seen only in relatively 
nearby samples is not stable which reflects strong fluctuations in the mean density. This 
is clearly seen in panel (a) where we show ^(s) for galaxies in the luminosity range 
— 19.5 < M < —18.5 in volumes ranging from 38 to 60 h^^ Mpc in depth as indicated in 
the figure. For the faintest range considered the volumes are relatively small and one finds 
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some variation in the correlation function at large separations with the amplitude of ^(s) 
computed in the 60 Mpc volume being significantly larger for s > 6 h"^ Mpc . This is 
not surprising since the size of the samples are comparable to the scale of inhomogeneities 
leading to important fluctuations in the mean density. This is not the case for the brighter 
subsamples probing volumes with depths D > 60 Mpc . In fact, although one can see 
in figure 1 (a) an increase in the amplitude of the correlation function between, for instance, 
samples D38 and D60, no significant luminosity segregation is detected when comparing the 
^(s) for different luminosity classes within the same volume. Instead ^(s) is seen to vary 
significantly with the depth of the volume considered. 

By contrast, for galaxies brighter than M — —19.5, shown in panels (b) and (c) there 
is no evidence for variations in the correlation strength with the depth of the sample. 
This suggests that sampling variations affect our results only for samples that probe small 
volumes of space. This conclusion is confirmed by the fact that the ratio of ^(s) for galaxies 
in the D91 subsample and for galaxies in the sample D138 is consistent with the ratio of 
^(s) for the same luminosity classes within the same volume. This indicates that luminosity 
segregation is responsible for the increase in the amplitude of ^ (s) between sample D91 
and D138 seen in Figure 1. We assume that the same is valid for galaxies brighter than 
M = —21, although we cannot verify this claim because the number of such bright galaxies 
rapidly decreases for smaller volumes. Nevertheless, 

we have checked that the mean density within this sample, calculated in shells of 
different depths, does not show any systematical variation. In particular, no drop of the 
density with depth has been detected. Therefore the clustering properties of the very bright 
galaxies previously shown are not a spurious effect due to an hypothetical incompletness. 
The properties of the subsample D168 are examined in detail in a forthcoming paper (Cappi 
et al.l996). 
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Based on the previous discussion we conclude that the variation of the amphtude of the 
correlation function observed in Figure 1 has two contributions. The first due to sampling 
fluctuations which affects primarily the clustering of faint galaxies probing small volumes. 
This affects samples roughly to D74 or galaxies fainter than M — —19. The second effect is 
consistent with the interpretation that the clustering of galaxies depends on luminosity. In 
the discussion below, we use the correlation function of volume-limited samples including all 
galaxies brighter than a limiting value to derive the relative bias of galaxies as a function of 
luminosity. We do that because, although useful to investigate the onset of the segregation, 
comparison of galaxies in magnitude bins are not suitable for comparison with theoretical 
predictions. 



3.3. Redshift Distortion and Morphology Effects 

As different structures are probed by the samples used in figure 1, variations in the 
amphtude of the correlation could also be induced by redshift distortions. On small scales, 
virialized systems tend to lower the correlation amplitude, while on intermediate scales 
streaming motions enhance the correlation amplitude, when measured in redshift space 
(Kaiser 1987). However, as da Costa et al. (1995) have shown that redshift distortions are 
relatively weak within the volume surveyed by the SSRS2, we do not expect a strong impact 
on the statistics. Nevertheless, in order to quantify this effect, we have checked the behavior 
of the correlation function in real space on the previously tested samples. We have verified 
that the trend with luminosity is the same, and that the amplification factor between the 
various samples is quantitatively the same in real space as in redshift space. This shows 
that redshift space distortion are not responsible for the luminosity biasing evidenced by 
SSRS2 data. 

We note that our results cannot be due to a morphology-luminosity relation since 
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Marzke et al. (1994) and da Costa et al. (1995) have shown that the luminosity functions for 
different morphological types are similar, except for irregulars. Moreover, we have verified 
that the morphological composition of our volume-limited samples is independent of the 
luminosity class adopted (Table 1, column (11)). 



3.4. Compcirison with previous results 

The above results give further support to conclusions reached by Hamilton (1988), 
Davis et al. (1988) and Lachieze-Rey, da Costa and Maurogordato (1992) from their 
analyses of the CfAl and SSRS redshift surveys, with significantly better statistics and 
extending the range of magnitudes over which the effect can be studied. First, both 
sampling effects and stronger clustering of more luminous galaxies are responsible for the 
variation of the correlation length with the depth of the sample. Sampling effects arc 
important for the smaller volumes and are difficult to disentangle from luminosity effects 
below M ~ —19. However, luminosity segregation at the low-luminosity end, if present at 
all, should be small, as figure 2 shows no evidence for a variation of the clustering amplitude 
of galaxies in the magnitude interval —18.5 < M < —17.5 relative to those in the interval 
—19.5 < M < —18.5 within the volume of depth SOh~^ Mpc . Second, for brighter galaxies, 
probing larger volumes, sampling effects are not important and for a given luminosity limit 
there is no significant variation of the correlation amplitude with depth. On the other hand, 
comparison of different luminosity classes within the same volume shows robust evidence 
for a relative increase in the clustering strength as a function of the luminosity. The relative 
amplitude is the same as that deduced from considering galaxies in the magnitude bins 
in the same volume. Third, the values of the relative bias when comparing the brightest 
galaxies in the sample to galaxies fainter that L* are in good agreement with the values 
obtained by Hamilton (1988) from the CfAl. 
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Signatures of luminosity segregation have also been detected by Park et al. (1994) 
from their power spectrum analysis of volume-limited subsamples drawn from CfA2. In 
particular, the amplification factor of ~ 1.4 between their subsamples at 101/i~^ Mpc, 
M < —19.7 and 130h~^ Mpc, M < —20.3 is consistent with our findings. However, their 
comparison of the power spectrum from bright/faint galaxies in volume- limited samples is 
not easily compared to our results as they use different magnitude bins. 

Evidence for luminosity-dependence of clustering have also been detected by Loveday 
et al. (1995) from a recent analysis of the correlation function on the Stromlo— APM redshift 
survey. Their conclusions differ from our findings as they show that while galaxies brighter 
than (—20 < M < —19) are more clustered than fainter ones (—19 < M < —15), there is 
no evidence for stronger clustering of galaxies with L > (—22 < M < —20). The reasons 
for these differences are at the present unclear. 



4. Bias and Compcirison to Models 

In order to quantify the dependence of the clustering strength on luminosity we define 
the relative bias b/b* of galaxies brighter than M towards galaxies brighter than M* 
(~ —19.5) as a function of scale s: 



We have adopted the following notation for the various statistics apphed: no subscript 
refers to galaxies brighter than M, and * subscript refers to galaxies brighter than M*. 

On those scales where the correlation function is well fitted by a power law, Js{s)/V 
and cr^(s) are proportional to ^(s), and therefore the ratio b/b^{s) can alternatively be 
expressed in terms of ^3(5) or of cr^(s) as 
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Ms) 



(10) 




(11) 



It is more convenient to work with expression (10) and (11) rather than with expression 
(9) as ^3(5) or cr^(s), which are integrated statistics, make possible the derivation of 
smoother ratios. 

As mentioned above, while the use of magnitude bins is useful to detect and to 
determine the onset of luminosity segregation, for comparison with theoretical predictions 
it is more convenient to use volume and absolute magnitude limited samples like those of 
Table 1. Results relative to M > —19.0 should be considered with caution before further 
interpretation, as sampling effects have been shown to be important for the smaller samples. 



The ratio 6/6* as a function of scale is presented in figure 4 as derived from ratios of 
73(5) (panel a) and o"^(s) (panel b). It is worth noting that the ratio is roughly constant in 
the interval 2 < s < 13h^^ Mpc for the subsamples considered, indicating that the relative 
bias is independent of scale, which makes it possible to define a mean linear biasing factor, 
which varies from 0.8 for the faintest sample (M < —19) to 2 for the brightest sample 
(M < -21.0). 

The variation of the linear bias with luminosity is shown in figure 5, where we show for 
comparison the same curve obtained by Hamilton (1988) from the CfAl survey. As can be 
seen there is a remarkable agreement at the bright end between these estimates based on 
samples probing different structures. It is important to emphasize that the SSRS2 sample 



EDITOR: PLACE FIGURE | HERE. 
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is not affected by uncertainties in the Virgo infall which has been used as an argument 
against Hamilton's interpretation of his results based on CfAl data. 

EDITOR: PLACE FIGURE | HERE. 

Luminosity biasing is expected to be present in the galaxy distribution according to 
most scenarios of galaxy formation. In the standard scheme of biased galaxy formation 
(Bardeen et al. 1986), galaxies are assumed to form from peaks above some global threshold 
in the primordial smoothed density field . This implies a dependence of the correlation 
function on the threshold which can be converted into a dependence on luminosity once 
a luminosity or circular velocity-peak relation is determined. Direct evidence that the 
amplitude of the correlation function increases with the circular velocity of the halos has 
been found both in numerical simulations of hierarchical clustering (White et al 1987) 
and in an extension of the Press-Schechter formalism (Mo and White 1995). Comparing 
these predictions to our data requires the knowledge of the relation between halos and 
galaxies. We have compared our results to these predictions by expressing circular velocities 
in terms of luminosities by the Tully-Fisher and Faber- Jackson relations (White, Tully 
and Davis 1988 and Valls-Gabaud, Alimi and Blanchard 1989). We have normalized the 
results for all the samples tested to the sample with galaxies brighter than L*. As shown 
in figure 5, even in the case of low normalization (cr8=0.5), it is difficult to mimic the 
steep increase of clustering present in the data at L > L*. Moreover, all models exhibit 
a clear luminosity dependence also at low luminosities which is difficult to test from the 
data as low-luminosity samples are also nearby and have been shown to be sensitive to 
sampling fluctuations. Still, as discussed before, the lack of luminosity segregation inside 
single volume-limited nearby samples as shown in figure 3 tends to indicate that a strong 
segregation at the low-luminosity end is quite unlikely. Note that this comparison is only 
suggestive as a one-to-one correspondence has been assigned to galaxy-halo pairs. More 
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elaborated modeling of the properties of galaxies and halos (Kauffmann, 1995) will allow a 
more precise comparison of luminosity biasing to the prediction of models for galaxies. 

An alternative view on the origin of biasing can be obtained from a description of the 
density field behavior at scales where it is strongly non-linear (Schaeffer 1984, 1985, Balian 
& Schaeffer 1989). This approach yields an analytic description of the matter fluctuations, 
based on an assumption about the scaling of the many-body matter correlation functions 
inspired by dynamical considerations (Peebles 1980, Balian & Schaeffer 1989). These models 
imply a bias in the distribution of high density spots. However, this bias is not acquired at 
the birth of the over-densities as a result of the initial Gaussian conditions, and then frozen. 
It is the result (Schaeffer 1985, 1987) of the non-hnear evolution, and is determined once the 
density fluctuations are known. A specific calculation (Bernardeau & Schaeffer 1992) shows 
that for all models with the above scaling, the bias depends only on the internal properties 
of the objects whose correlations are considered. More specifically it increases as a function 
of a unique parameter, closely related to the gravitational potential, and therefore increases 
with luminosity. This analytical model gives explicitly the asymptotic behavior for the 
brightest galaxies: h/b* is nearly proportional to L/L*. It is very successful in reproducing 
the steep increase of the correlation function for bright galaxies L > L* but it predicts a 
strong luminosity segregation at faint luminosities which is very difficult to reconcile with 
the data (figure 5). Larger samples of galaxies going to fainter magnitudes are therefore 
necessary to test the models at the low-luminosity end. 

5. Summciry 

Using the recently completely SSRS2 redshift survey we have investigated the 
dependence of the clustering properties of galaxies on luminosity. The SSRS2 sample is 
particularly suitable for this study because of the large number of galaxies, the paucity of 
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rich clusters within the volume and the small amplitude of redshift distortions induced by 
peculiar motions both at small and large scales. Analysis of volume-limited samples using 
both the two-point correlation function and the variance of counts in cells indicate that 
galaxy clustering depends on luminosity. The effect is small for faint galaxies but increases 
dramatically for galaxies brighter than L*. The correlation length varies from about 5 
Mpc , for the faintest galaxies up to 8 h^^ Mpc for the galaxies brighter than M = —20.5. 
For the brightest galaxies in the sample with M < —21 the correlation length reaches values 
as high as ~ 16 Mpc . These galaxies are not preferentially found in rich clusters and 
the study of their properties will be the subject of a separate paper. The relative bias is 
scale independent up to ~ 13 Mpc and the linear bias between the brightest and 
galaxies is in the range 1.5-2. 

A preliminary comparison with presently available models of galaxy formation shows 
serious difficulties in modeling simultaneously the faint and bright ends of the bias variation 
with luminosity evidenced in the SSRS2. However, more data is necessary to allow a better 
determination of the behavior of luminosity biasing at faint luminosities which is strongly 
affected by sampling variations for nearby catalogs, although it is constrained to be small 
by our analysis. On the theoretical side, a better understanding galaxy formation and of the 
galaxy-dark-halo connection is required in order to use these results as useful constraints to 
galaxy formation models. 

We would like to thank all the SSRS2 collaborators for allowing us to use the data in 
advance of its publication. We are grateful to George Efstathiou, Guinevere Kauffmann, 
Amos Yahil and Simon White for fruitful discussions. LNdC would like to thank the 
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Fig. 1. — The two-point correlation function in redshift space ^(s) (panel a), the Jsis) 
integral (panel b) and the variance o'^(s) (panel c) calculated from counts in cells for SSRS2 
volume-limited subsamples identified as follows : D30 (stars); D38 (full triangles); D48 
(opened triangles): D60 (full pentagons): D74 (opened pentagons); D91 (full squares); D112 
(open squares); D138 (full circles); D168 (open circles). 

Fig. 2. — The variation of ^(s) with luminosity is tested within subsamples at a given depth. 
In each panel, wc present ^(s) for two luminosity classes of one magnitude range within the 
same volume. In panel a), galaxies with —18.5 < M < —17.5 and —19.5 < M < —18.5 are 
compared within a volume of 3Sh~^Mpc depth. In panel b), galaxies with —19.5 < M < 
—18.5 and —20.5 < M < —19.5 are compared within a volume of 60h~^Mpc depth. In panel 
c), galaxies with —20.5 < M < —19.5 and —21.5 < M < —20.5 are compared within a 
volume of 91h~^Mpc depth. 

Fig. 3. — The redshift space correlation function ^(s) is displayed in various panels, each 
corresponding to a given luminosity class: —19.5 < M < —18.5 in panel a), —20.5 < M < 
—19.5. in panel b), —21.5 < M < —20.5 in panel c). Within each panel, ^(s) is presented 
for three different volumes. Symbols vary from opened triangles to opened squares and full 
squares with increasing size of the sample. 

Fig. 4. — The biasing factor as a function of separation, for subsamples D74 to D168. It is 
normalized to the value b* calculated from sample D91 corresponding to galaxies brighter 
than L*. In panel a), the biasing factor is derived from J3, and in panel b) from a^. 

Fig. 5. — Evolution of b/b* (averaged between 2h~^Mpc and ISh'^Mpc) with the limiting 
absolute magnitude of the sample . Our results from SSRS2 are reported as full squares 
when derived from J3, and as open squares when derived from a^. CfA results of Hamilton 
(1988) are shown as opened triangles. Results from CDM simulations of White et al. (1987) 
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as given by Valls-Gabaud et al. (1989) are plotted as a long-dashed line. On this graph are 
over-plotted analytical predictions by Mo and White (1995), with low (o-g = 0.5, dotted hne) 
and high normalizations (erg = 1, short-dashedline) and by Bernardeau & Schaeffer (1992) 
(dashed-dotted line). 
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Table 1. Subsamples of the SSRS2. 



Sample 


L) 






n 


so 


7 


so' 


Interval 


_2 


TIT? 




(/i-^Mpc) 






(10-3/l3Mpc-3) 


(ft-iMpc) 




(fe-iMpc) 


(/i-iMpc) 






D30 


30.69 


-17.0 


301 


27.2 


4.3±0.5 


1.7±0.2 


4.3±0.4 


2.0-5.2 




2.3 


D38 


38.35 


-17.5 


329 


15.5 


4.5±0.4 


1.6±0.2 


4.4±0.3 


2.0-5.2 




2.2 


D48 


47.83 


-18.0 


487 


11.8 


3.5±0.2 


2.1±0.2 


3.5±0.2 


2.0-5.0 


0.3 


2.2 


D60 


59.52 


-18.5 


774 


9.7 


4.8±0.5 


1.3±0.1 


4.6±0.3 


2.0-8.0 


0.5 


2.0 


D74 


73.86 


-19.0 


770 


5.1 


5.5±0.4 


1.4±0.1 


5.2±0.4 


2.0-8.0 


0.7 


1.9 


D91 


91.37 


-19.5 


755 


2.6 


5.8±0.6 


1.7±0.2 


5.8±0.4 


2.0-12.6 


0.9 


2.1 


D112 


112.60 


-20.0 


593 


1.1 


5.5±0.6 


2.0±0.3 


5.5±0.4 


2.0-12.6 


0.9 


2.2 


D138 


138.12 


-20.5 


268 


0.3 


8.0±0.4 


2.0±0.2 


8.5±0.5 


2.0-15.8 


1.5 


2.5 


D168 


168.52 


-21.0 


133 


0.1 


15.8±2.9 


1.5±0.3 


15.2±2.0 


2.5-25.0 


3.0 


1.7 



Table 2. Characteristics of subsamples in magnitude bins. 



Sample 


D 






(/i-iMpc) 




-18.5/-17.5 


38.35 


187 


-19.5/-18.5 


38.35 


105 


-19.5/-18.5 


47.83 


227 


-19.5/-18.5 


59.52 


563 


-20.5/-19.5 


59.52 


193 


-20.5/-19.5 


73.86 


329 


-20.5/-19.5 


91.37 


687 


-21.5/-20.5 


91.37 


67 


-21.5/-20.5 


112.60 


155 


-21.5/-20.5 


138.12 


268 
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Fig. 5. — 



